Cell-type-specific mRNA transcription and degradation kinetics in zebrafish embryogenesis from metabolically labeled single-cell RNA-seq.

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the relative contributions of mRNA transcription and degradation to shaping those profiles remains challenging, especially within embryos with diverse cellular identities. Here, we combine single-cell RNA-Seq and metabolic labeling to capture temporal cellular transcriptomes of zebrafish embryos where newly-transcribed (zygotic) and pre-existing (maternal) mRNA can be distinguished. We introduce kinetic models to quantify mRNA transcription and degradation rates within individual cell types during their specification. These models reveal highly varied regulatory rates across thousands of genes, coordinated transcription and destruction rates for many transcripts, and link differences in degradation to specific sequence elements. They also identify cell-type-specific differences in degradation, namely selective retention of maternal transcripts within primordial germ cells and enveloping layer cells, two of the earliest specified cell types. Our study provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development.

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Germ cells do not progress through spermatogenesis in the infertile zebrafish testis.

Vertebrate spermatogonial stem cells maintain sperm production over the lifetime of an animal but fertility declines with age. While morphological studies have greatly informed our understanding of typical spermatogenesis, the molecular and cellular mechanisms underlying spermatogenesis are not yet understood, particularly with respect to the onset of fertility. We used single-cell RNA sequencing to generate a developmental atlas of the zebrafish testis. Using 5 timepoints across the adult life of a zebrafish, we described cellular profiles in the testis during and after fertility. While all germ cell stages of spermatogenesis are detected in testes from fertile adult zebrafish, testes from older infertile males only contained spermatogonia and a reduced population of spermatocytes. These remaining germ cells are transcriptionally distinct from fertile spermatogonia. Immune cells including macrophages and lymphocytes drastically increase in abundance in infertile testes. Our developmental atlas reveals the cellular changes as the testis ages and defines a molecular roadmap for the regulation of male fertility.

Single-cell temporal dynamics reveals the relative contributions of transcription and degradation to cell-type specific gene expression in zebrafish embryos.

During embryonic development, pluripotent cells assume specialized identities by adopting particular gene expression profiles. However, systematically dissecting the underlying regulation of mRNA transcription and degradation remains a challenge, especially within whole embryos with diverse cellular identities. Here, we collect temporal cellular transcriptomes of zebrafish embryos, and decompose them into their newly-transcribed (zygotic) and pre-existing (maternal) mRNA components by combining single-cell RNA-Seq and metabolic labeling. We introduce kinetic models capable of quantifying regulatory rates of mRNA transcription and degradation within individual cell types during their specification. These reveal different regulatory rates between thousands of genes, and sometimes between cell types, that shape spatio-temporal expression patterns. Transcription drives most cell-type restricted gene expression. However, selective retention of maternal transcripts helps to define the gene expression profiles of germ cells and enveloping layer cells, two of the earliest specified cell-types. Coordination between transcription and degradation restricts expression of maternal-zygotic genes to specific cell types or times, and allows the emergence of spatio-temporal patterns when overall mRNA levels are held relatively constant. Sequence-based analysis links differences in degradation to specific sequence motifs. Our study reveals mRNA transcription and degradation events that control embryonic gene expression, and provides a quantitative approach to study mRNA regulation during a dynamic spatio-temporal response.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development.

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 hours post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and suggest new long-term cycling populations. Focused analyses of non-skeletal muscle and the endoderm identified transcriptional profiles of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and homologs of recently discovered human best4+ enterocytes. The transcriptional regulators of these populations remain unknown, so we reconstructed gene expression trajectories to suggest candidates. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Differentiation trajectories of the Hydra nervous system reveal transcriptional regulators of neuronal fate.

The small freshwater cnidarian polyp Hydra vulgaris uses adult stem cells (interstitial stem cells) to continually replace neurons throughout its life. This feature, combined with the ability to image the entire nervous system (Badhiwala et al., 2021; Dupre & Yuste, 2017) and availability of gene knockdown techniques (Juliano, Reich, et al., 2014; Lohmann et al., 1999; Vogg et al., 2022), makes Hydra a tractable model for studying nervous system development and regeneration at the whole-organism level. In this study, we use single-cell RNA sequencing and trajectory inference to provide a comprehensive molecular description of the adult nervous system. This includes the most detailed transcriptional characterization of the adult Hydra nervous system to date. We identified eleven unique neuron subtypes together with the transcriptional changes that occur as the interstitial stem cells differentiate into each subtype. Towards the goal of building gene regulatory networks to describe Hydra neuron differentiation, we identified 48 transcription factors expressed specifically in the Hydra nervous system, including many that are conserved regulators of neurogenesis in bilaterians. We also performed ATAC-seq on sorted neurons to uncover previously unidentified putative regulatory regions near neuron-specific genes. Finally, we provide evidence to support the existence of transdifferentiation between mature neuron subtypes and we identify previously unknown transition states in these pathways. All together, we provide a comprehensive transcriptional description of an entire adult nervous system, including differentiation and transdifferentiation pathways, which provides a significant advance towards understanding mechanisms that underlie nervous system regeneration.

Gene module reconstruction elucidates cellular differentiation processes and the regulatory logic of specialized secretion.

During differentiation, cells become structurally and functionally specialized, but comprehensive views of the underlying remodeling processes are elusive. Here, we leverage scRNA-seq developmental trajectories to reconstruct differentiation using two secretory tissues as a model system - the zebrafish notochord and hatching gland. First, we present an approach to integrate expression and functional similarities for gene module identification, revealing dozens of gene modules representing known and newly associated differentiation processes and their temporal ordering. Second, we focused on the unfolded protein response (UPR) transducer module to study how general versus cell-type specific secretory functions are regulated. By profiling loss- and gain-of-function embryos, we found that the UPR transcription factors creb3l1, creb3l2, and xbp1 are master regulators of a general secretion program. creb3l1/creb3l2 additionally activate an extracellular matrix secretion program, while xbp1 partners with bhlha15 to activate a gland-specific secretion program. Our study offers a multi-source integrated approach for functional gene module identification and illustrates how transcription factors confer general and specialized cellular functions.

Emergence of Neuronal Diversity during Vertebrate Brain Development.

Neurogenesis comprises many highly regulated processes including proliferation, differentiation, and maturation. However, the transcriptional landscapes underlying brain development are poorly characterized. We describe a developmental single-cell catalog of ∼220,000 zebrafish brain cells encompassing 12 stages from embryo to larva. We characterize known and novel gene markers for ∼800 clusters and provide an overview of the diversification of neurons and progenitors across these time points. We also introduce an optimized GESTALT lineage recorder that enables higher expression and recovery of Cas9-edited barcodes to query lineage segregation. Cell type characterization indicates that most embryonic neural progenitor states are transitory and transcriptionally distinct from neural progenitors of post-embryonic stages. Reconstruction of cell specification trajectories reveals that late-stage retinal neural progenitors transcriptionally overlap cell states observed in the embryo. The zebrafish brain development atlas provides a resource to define and manipulate specific subsets of neurons and to uncover the molecular mechanisms underlying vertebrate neurogenesis.

Stem cell differentiation trajectories in Hydra resolved at single-cell resolution.

The adult Hydra polyp continually renews all of its cells using three separate stem cell populations, but the genetic pathways enabling this homeostatic tissue maintenance are not well understood. We sequenced 24,985 Hydra single-cell transcriptomes and identified the molecular signatures of a broad spectrum of cell states, from stem cells to terminally differentiated cells. We constructed differentiation trajectories for each cell lineage and identified gene modules and putative regulators expressed along these trajectories, thus creating a comprehensive molecular map of all developmental lineages in the adult animal. In addition, we built a gene expression map of the Hydra nervous system. Our work constitutes a resource for addressing questions regarding the evolution of metazoan developmental processes and nervous system function.

Single-cell reconstruction of developmental trajectories during zebrafish embryogenesis.

During embryogenesis, cells acquire distinct fates by transitioning through transcriptional states. To uncover these transcriptional trajectories during zebrafish embryogenesis, we sequenced 38,731 cells and developed URD, a simulated diffusion-based computational reconstruction method. URD identified the trajectories of 25 cell types through early somitogenesis, gene expression along them, and their spatial origin in the blastula. Analysis of Nodal signaling mutants revealed that their transcriptomes were canalized into a subset of wild-type transcriptional trajectories. Some wild-type developmental branch points contained cells that express genes characteristic of multiple fates. These cells appeared to trans-specify from one fate to another. These findings reconstruct the transcriptional trajectories of a vertebrate embryo, highlight the concurrent canalization and plasticity of embryonic specification, and provide a framework with which to reconstruct complex developmental trees from single-cell transcriptomes.

Spatial reconstruction of single-cell gene expression data.

Spatial localization is a key determinant of cellular fate and behavior, but methods for spatially resolved, transcriptome-wide gene expression profiling across complex tissues are lacking. RNA staining methods assay only a small number of transcripts, whereas single-cell RNA-seq, which measures global gene expression, separates cells from their native spatial context. Here we present Seurat, a computational strategy to infer cellular localization by integrating single-cell RNA-seq data with in situ RNA patterns. We applied Seurat to spatially map 851 single cells from dissociated zebrafish (Danio rerio) embryos and generated a transcriptome-wide map of spatial patterning. We confirmed Seurat's accuracy using several experimental approaches, then used the strategy to identify a set of archetypal expression patterns and spatial markers. Seurat correctly localizes rare subpopulations, accurately mapping both spatially restricted and scattered groups. Seurat will be applicable to mapping cellular localization within complex patterned tissues in diverse systems.

Response to Nodal morphogen gradient is determined by the kinetics of target gene induction.

Morphogen gradients expose cells to different signal concentrations and induce target genes with different ranges of expression. To determine how the Nodal morphogen gradient induces distinct gene expression patterns during zebrafish embryogenesis, we measured the activation dynamics of the signal transducer Smad2 and the expression kinetics of long- and short-range target genes. We found that threshold models based on ligand concentration are insufficient to predict the response of target genes. Instead, morphogen interpretation is shaped by the kinetics of target gene induction: the higher the rate of transcription and the earlier the onset of induction, the greater the spatial range of expression. Thus, the timing and magnitude of target gene expression can be used to modulate the range of expression and diversify the response to morphogen gradients.

From egg to gastrula: how the cell cycle is remodeled during the Drosophila mid-blastula transition.

Many, if not most, embryos begin development with extremely short cell cycles that exhibit unusually rapid DNA replication and no gap phases. The commitment to the cell cycle in the early embryo appears to preclude many other cellular processes that only emerge as the cell cycle slows just prior to gastrulation at a major embryonic transition known as the mid-blastula transition (MBT). As reviewed here, genetic and molecular studies in Drosophila have identified changes that extend S phase and introduce a post-replicative gap phase, G2, to slow the cell cycle. Although many mysteries remain about the upstream regulators of these changes, we review the core mechanisms of the change in cell cycle regulation and discuss advances in our understanding of how these might be timed and triggered. Finally, we consider how the elements of this program may be conserved or changed in other organisms.

Mechanism and regulation of Cdc25/Twine protein destruction in embryonic cell-cycle remodeling.

BACKGROUND: In Drosophila embryos, the midblastula transition (MBT) dramatically remodels the cell cycle during the 14(th) interphase. Before the MBT, each cycle is composed of only a short S phase and mitosis. At the MBT, S phase is dramatically lengthened by the onset of late replication, and a G2 phase is introduced. Both changes set the stage for gastrulation and require downregulation of Cdc25 phosphatase, which was previously attributed to the elimination of its transcripts at the MBT. RESULTS: Premature removal of cdc25 transcripts by RNAi did not affect progression to the MBT. Instead, an antibody against the Cdc25 isoform Twine showed that Twine protein was abundant and stable until the MBT, when it was destabilized and rapidly eliminated. Persistence of pre-MBT levels of Twine was sufficient to prevent cell-cycle slowing. Twine protein destruction was timed by the nucleocytoplasmic ratio and depended on the activation of zygotic transcription at the MBT, including expression of the gene tribbles, whose activity was sufficient to trigger Twine destruction and was required for prompt Twine disappearance. CONCLUSIONS: We propose that the developmentally regulated destruction of Twine protein is a critical switch that contributes to the cell-cycle change at the MBT, including the addition of a G2 phase and onset of late replication. Moreover, we show that this destruction is triggered by the nucleocytoplasmic ratio-dependent onset of zygotic transcription of tribbles and other unknown genes.

Different cyclin types collaborate to reverse the S-phase checkpoint and permit prompt mitosis.

Precise timing coordinates cell proliferation with embryonic morphogenesis. As Drosophila melanogaster embryos approach cell cycle 14 and the midblastula transition, rapid embryonic cell cycles slow because S phase lengthens, which delays mitosis via the S-phase checkpoint. We probed the contributions of each of the three mitotic cyclins to this timing of interphase duration. Each pairwise RNA interference knockdown of two cyclins lengthened interphase 13 by introducing a G2 phase of a distinct duration. In contrast, pairwise cyclin knockdowns failed to introduce a G2 in embryos that lacked an S-phase checkpoint. Thus, the single remaining cyclin is sufficient to induce early mitotic entry, but reversal of the S-phase checkpoint is compromised by pairwise cyclin knockdown. Manipulating cyclin levels revealed that the diversity of cyclin types rather than cyclin level influenced checkpoint reversal. We conclude that different cyclin types have distinct abilities to reverse the checkpoint but that they collaborate to do so rapidly.

Embryonic onset of late replication requires Cdc25 down-regulation.

The Drosophila midblastula transition (MBT), a major event in embryogenesis, remodels and slows the cell cycle. In the pre-MBT cycles, all genomic regions replicate simultaneously in rapid S phases that alternate with mitosis, skipping gap phases. At the MBT, down-regulation of Cdc25 phosphatase and the resulting inhibitory phosphorylation of the mitotic kinase Cdk1 create a G2 pause in interphase 14. However, an earlier change in interphase 14 is the prolongation of S phase. While the signals modifying S phase are unknown, the onset of late replication-where replication of constitutively heterochromatic satellite sequences is delayed-extends S-phase 14. We injected Cdc25 mRNA to bypass the developmentally programmed down-regulation of Cdc25 at the MBT. Introduction of either Cdc25 isoform (String or Twine) or enhanced Cdk1 activity triggered premature replication of late-replicating sequences, even after their specification, and thereby shortened S phase. Reciprocally, reduction of Cdk1 activity by knockdown of mitotic cyclins extended pre-MBT S phase. These findings suggest that high Cdc25 and Cdk1 contribute to the speed of the rapid, pre-MBT S phases and that down-regulation of these activities plays a broader role in MBT-associated changes than was previously suspected.

Influence of cyclin type and dose on mitotic entry and progression in the early Drosophila embryo.

Cyclins are key cell cycle regulators, yet few analyses test their role in timing the events that they regulate. We used RNA interference and real-time visualization in embryos to define the events regulated by each of the three mitotic cyclins of Drosophila melanogaster, CycA, CycB, and CycB3. Each individual and pairwise knockdown results in distinct mitotic phenotypes. For example, mitosis without metaphase occurs upon knockdown of CycA and CycB. To separate the role of cyclin levels from the influences of cyclin type, we knocked down two cyclins and reduced the gene dose of the one remaining cyclin. This reduction did not prolong interphase but instead interrupted mitotic progression. Mitotic prophase chromosomes formed, centrosomes divided, and nuclei exited mitosis without executing later events. This prompt but curtailed mitosis shows that accumulation of cyclin function does not directly time mitotic entry in these early embryonic cycles and that cyclin function can be sufficient for some mitotic events although inadequate for others.